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The nonlinear Breit-Wheeler process is studied in the presence of strong and short laser pulses. 

We show that for a relativistically intense plane-wave laser field many features of the momentum dis¬ 
tribution of the produced electron-positron pair like its extension, region of highest probability and 
carrier-envelope phase effects can be explained from the classical evolution of the created particles 
in the background field. To this end an intuitive semiclassical picture based on the local constant- 
crossed field approximation applied on the probability-amplitude level is established and compared 
with the standard approach used in QED-PIC codes. The main difference is the substructure of 
the spectrum, which results from interference effects between macroscopically separated formation 
regions. In order to compare the predictions of the semiclassical approach with exact calculations, a 
very fast numerical scheme is introduced. It renders the calculation of the fully differential spectrum 
on a grid which resolves all interference fringes feasible. Finally, the difference between classical and 
quantum absorption of laser four-momentum in the process is pointed out and the dominance of the 
former is proven. As a self-consistent treatment of the quantum absorption is not feasible within 
existing QED-PIC approaches, our results provide reliable error estimates relevant for regimes where 
the laser depletion due to a developing QED cascade becomes significant. 


I. INTRODUCTION 

From a conceptual point of view the transformation of 
light into matter is one of the most appealing physical 
processes. The possibility to create an electron-positron 
pair by merging two real photons (Breit-Wheeler process 
[1]) is a direct manifestation of the equivalence of mass 
and energy, postulated first by Einstein. Hitherto, how¬ 
ever, this process has not been observed in a laboratory. 
Experimentally, electron-positron photoproduction was 
studied only indirectly via the trident process, i.e. by col¬ 
liding a highly relativistic electron beam with an optical 
laser in the SLAC E-144 experiment [2] (the experimental 
findings could be explained within the two-step approx¬ 
imation, i.e. by the assumption that the electrons first 
radiate real gamma photons which subsequently decay 
into pairs via the Breit-Wheeler process). To understand 
the results of the measurement theoretically, nonlinear 
effects (i.e. the simultaneous absorption of several laser 
photons) must be taken into account [2-5]. 

Due to the continuous improvement of laser technol¬ 
ogy the experimental observation of the nonlinear gener¬ 
alization of the Breit-Wheeler process (see Fig. 1) is now 
within reach, e.g., at upcoming high-power laser facili¬ 
ties like the APOLLON-lOP laser [9], the Extreme Light 
Infrastructure [10], the Vulcan 10PW laser [11] at the 
Central Laser Facility [12], and the Exawatt Center for 
Extreme Light Studies [13]. Therefore, this process has 
recently been considered by several authors [14-30] (see 
also the reviews [31-38]). 

The decay of a photon into an electron-positron pair 
is an intrinsic quantum process, which has no classical 
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analogue and must be described in the realm of quan¬ 
tum field theory, e.g., by calculating the corresponding 
S'-matrix element. This implies that we can only deter¬ 
mine the total probability for the decay and the asymp¬ 
totic momentum distribution for the final particles. 

However, it is well known that in the case of a plane- 
wave laser field with electric field amplitude E 0 and cen¬ 
tral angular frequency to the formation region of the basic 
QED processes nonlinear Compton scattering and non¬ 
linear Breit-Wheeler pair production are £-times smaller 
than the laser period in the ultra-relativistic regime 
£ 1. Here, £ = |e| E 0 /(mwc) is the classical inten¬ 

sity parameter, with e < 0 and m denoting the electron 
charge and mass, respectively [32, 37]. Hence, the total 
probability for nonlinear Compton scattering and non¬ 
linear Breit-Wheeler pair production can be calculated 
by applying the local constant-crossed field approxima¬ 
tion, i.e. by averaging the corresponding probability in a 



FIG. 1. Leading-order Feynman diagram for electron-positron 
photoproduction inside a plane-wave background field (non¬ 
linear Breit-Wheeler process). The double lines represent 
Volkov states (solutions of the interacting Dirac equation, 
which take the plane-wave background field into account ex¬ 
actly [6, 7]), the wiggly line the incoming photon. As long 
as the total pair-production probability is much smaller than 
unity, the spectrum for the final particles is determined to a 
good accuracy by simply evaluating this diagram and neglect¬ 
ing radiative corrections [8]. 
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constant-crossed field over the laser pulse [8, 32, 37, 39- 
43], 

As pointed out by Ritus [37], this procedure is justi¬ 
fied for the calculation of total probabilities but fails in 
general for the momentum distribution of the final parti¬ 
cles (this has also been recently observed numerically in 
[44] for nonlinear Compton scattering). The reason are 
interference effects arising from processes occurring at 
macroscopically separated space-time points, which are 
neglected from the beginning when the averaging proce¬ 
dure is applied to the probabilities. As a laser field is 
oscillatory, each cycle typically contains two formation 
regions where the electron-positron pair is created with 
the same asymptotic quantum numbers. An interference 
between different pathways is therefore expected from 
general principles similar to a multi-slit experiment. 

The fact that the substructure of the spectrum can be 
attributed to interference effects is well known from other 
processes. For nonlinear Thomson (Compton) scattering 
this was, e.g., reported in Refs. [45-49], for Schwinger 
pair production interference effects are, e.g., discussed in 
Refs. [50-53]. 

Nonetheless, the average over the laser pulse shape on 
the probability level (incoherent summation of all pos¬ 
sible processes) rather than the amplitude level (coher¬ 
ent summation of all possible processes) is the state-of- 
the-art approach for the implementation of strong-field 
QED processes in so-called particle-in-cell (PIC) codes 
[54-81]. Therefore, it is desirable to revisit this approach 
and show how it could be extended if a higher precision 
becomes necessary. To this end we study here for the first 
time the quantitative influence of interference effects on 
the spectra for electron-positron photoproduction in the 
semiclassical regime. 

By applying a stationary-phase analysis to the leading- 
order S'-matrix element for electron-positron photopro¬ 
duction it is shown below that in the strong-field regime 
£ 1 all qualitative features of the spectrum (includ¬ 

ing the substructure caused by interference effects) can 
be understood using the following semiclassical descrip¬ 
tion: at each laser phase the pair-production probability 
amplitude is calculated by employing the local constant- 
crossed held approximation. It predicts the local momen¬ 
tum spectrum of the pair immediately after the particles 
are brought on shell. The latter is then employed as 
initial condition for the classical propagation which pro¬ 
vides the asymptotic momentum distribution. Finally, 
the probability for the production of a pair with given 
asymptotic momenta is obtained by squaring the prob¬ 
ability amplitude, taking the interference between pairs 
which have the same asymptotic momenta but originat¬ 
ing from different formation regions into account. 

This intuitive picture provides a clear physical reason 
for many properties of the asymptotic momentum distri¬ 
bution like its extension and the regions of highest proba¬ 
bility. Moreover, it explains the strong dependence of the 
spectrum on the carrier-envelope phase (CEP) for ultra- 
short laser pulses reported in [22] . In the context of vac¬ 


uum pair production classical features have been identi¬ 
fied in [82] . Furthermore, a similar semiclassical analysis 
has been carried out in [46, 83, 84] to explain various 
features of the emission spectra for nonlinear Compton 
scattering and has been exploited in [84] to put forward 
a scheme for determining the CEP in ultra-short and 
ultra-intense laser beams via nonlinear single-Compton 
scattering. 

We point out that electron-positron photoproduction 
has a lot of commonalities with laser-induced ionization 
processes. In fact, the procedure outlined above is closely 
related to similar approaches used in atomic physics to 
describe the time evolution of an electron after tunnel 
ionization [85, 86]. 

Interestingly, this semiclassical three-step model allows 
us to distinguish between a “classical” and a “quantum” 
absorption of laser four-momentum in the process. As the 
decay of a photon into an electron-positron pair is forbid¬ 
den in vacuum, a certain amount of laser four-momentum 
must be absorbed initially to bring the massive particles 
on shell (this part will be called quantum absorption here, 
see Sec. HIE). Afterwards, the charged particles are fur¬ 
ther accelerated by the laser field, which implies a clas¬ 
sical energy-momentum transfer. Below, these two pro¬ 
cesses are distinguished for the first time and it is shown 
that classical absorption dominates for £ 1. Corre¬ 

spondingly, the laser is predominantly depleted from the 
classical energy-momentum transfer. 

To establish the validity of the outlined semiclassical 
approach, we compare its predictions with a full numeri¬ 
cal calculation of the leading-order S'-matrix element (see 
Fig. 1). In this way we show that already for £ > 5 the in¬ 
terference substructure obtained from the local constant- 
crossed field approximation applied on the probability 
amplitude level is in very good agreement with the full 
numerical result. 

As the S-matrix integrals are highly oscillating in the 
regime £ 1, the numerical calculation is, in princi¬ 

ple, a challenging task [21, 22]. We present here a new 
scheme, which is substantially more efficient than other 
employed methods (see Sec. IIC). Hence, it becomes fea¬ 
sible to evaluate the three-dimensional differential proba¬ 
bility on a grid which is fine enough to completely resolve 
the interference substructure of the spectrum. To obtain 
the total pair-production probability, however, the nu¬ 
merical evaluation of the compact expressions derived in 
[8] is still much faster. 

Note that for nonlinear Compton scattering the final 
phase space has been studied in [87-90] for moderate 
values of the parameter £ (£ < 10) and the computation 
of total probabilities can also be significantly simplified 
by employing the method proposed in [91]. 

The paper is organized as follows. In Sec. II our nota¬ 
tion is introduced and the well-known leading-order ex¬ 
pression for the pair-production probability [see Eq. (8)] 
is expressed using Lorentz-invariant momentum param¬ 
eters (see Sec. IIA for details). Furthermore, the nu¬ 
merical scheme for the numerical evaluation of the pair- 
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production probability is described in Sec. IIC. After¬ 
ward, the main new result of the present paper, the dis¬ 
entanglement between classical and quantum aspects of 
the pair-production process, is presented in Sec. III. 

From now on we use natural units K = c = 1 and 
Heaviside-Lorentz units for charge [a = e 2 /(47r) ss 1/137 
denotes the fine-structure constant], the notation agrees 
with [8, 92]. 


II. PAIR PRODUCTION PROBABILITY 

In this paper we consider the decay of a photon with 
four-momentum q^ ( q 2 = 0) and polarization four-vector 
into an electron and a positron with four-momentum 
Pi and p%, respectively (p 2 = m 2 ). In vacuum, this 
transition is forbidden by energy-momentum conserva¬ 
tion (for all photon energies). Inside an electromag¬ 
netic background field, however, the so-called (nonlinear) 
Breit-Wheeler process is in general allowed (see Fig. 1). 
Here, we focus on plane-wave laser fields described by the 
classical field tensor 

= E rmi /r = w 

i=1,2 


where cf> = kx is the laser phase, x^ the position four- 
vector, k M the characteristic four-momentum of the laser 
photons and a/ and (</>)] < 1 characterize the strength 
and the shape of the field, respectively, along the two 
possible polarization directions [i = 1,2; ^(± 00 ) = 
^(± 00 ) = 0; the prime denotes the derivative of a func¬ 
tion with respect to the argument]. Furthermore, we 
introduce the classical intensity parameters 

£, = sV-“?' t = P+$ C) 

and the quantum-nonlinearity parameter \ = ( kq/m 2 )l; 
(note that kq > 0 as q 2 = 0 here). 

If a laser field is sufficiently intense and depletion ef¬ 
fects are negligible, it can be considered as a classical 
background field. Correspondingly, the pair-production 
probability is obtainable using the formalism of strong- 
field QED in the Furry picture [93] (for more details see, 
e.g., [31-38, 94-96] and the literature cited in the intro¬ 
duction) . 


A. Invariant momentum parameters 

Using the canonical light-cone basis (see App. A for 
more details), we introduce the Lorentz-invariant mo¬ 
mentum parameters r, t\ and £2 


Pi = r'q M + a'k? + t'imAi + t^mA^, (3a) 

P 2 = — rq^ — tiviAi — t 2 «iA2. (3b) 


Here, = tj and r' = r + 1 due to momentum conserva¬ 
tion and the quantities 


s = 


1 m 2 
2r kq 


(1 + t\ + tl), 


1 to 2 
2 r' kq 


(1 + (4) 


are determined by the on-shell conditions p\ = p\ = m 2 . 
We point out that we consider only plane-wave fields with 
a finite duration. Correspondingly, the four-momentum 
q^ (pi, P 2 ) denotes the asymptotic four-momentum in 
vacuum before (after) the interaction with the laser pulse. 
In particular, we do not introduce dressed masses and 
momenta for the electron and the positron [32, 37]. 

For later convenience we also define the quantities [37] 


1 _ kpx - kp 2 _ kpi - kp 2 

2 2 kq 2 (kpi + fcp 2 ) ’ 

1 4 (fc?) 2 

W r(r + 1) 1 — 4i? 2 (kp\)(kp 2 ) 


(5a) 

(5b) 


to characterize how the initial photon four-momentum is 
split between the outgoing particles (unlike r and r' the 
parameters R and w are anti-symmetric and symmet¬ 
ric with respect to the electron and the positron four- 
momentum, respectively). As kq > 0, kpi > 0 and 
r' = r + 1 we conclude from Eq. (3) that r € (—1,0), 
r' € (0,1), R € (—1/2, +1/2) and w £ [4, 00 ). 

Note that the map from R -A w has no inverse (i.e. the 
information about the sign of R is lost). Correspondingly, 
the quantities kq and w [see Eq. (5)] specify kpi and fcp 2 
uniquely up to the sign of R. However, we will later see 
that this sign does not influence the spin-summed pair- 
production probability. 

Using the Lorentz-invariant momentum parameters r 
(R), t\ and f 2 to describe the asymptotic momenta of 
the created electron-positron pair has the advantage that 
they characterize the process in a frame-independent way. 
In particular, we do not have to work in a frame where 
the collision is head on. However, in this frame the pa¬ 
rameters ti have a simple interpretation, as they measure 
the transverse momentum of the pair in units of the elec¬ 
tron rest mass, whereas R measures in general how the 
initial photon four-momentum q^ is distributed between 
the electron and the positron. 

Note that in the seminal papers [37, 39] the differen¬ 
tial pair-production rate inside a constant-crossed field is 
expressed with respect to the two parameters u and r, 
which are related to the parameters w and t 2 introduced 
here as follows 


qf 2 q 


V(Plf 2 Pl)(.P2f 2 P2 ) 
Plf*P2 


= 4 M, 


tn = T 2 


m 




(6a) 

(6b) 


(/ M!y = fi")- Therefore, the comparison between the 
results obtained here for a laser pulse and those reported 
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previously for a constant-crossed field is straightforward 
in canonical light-cone coordinates. 

The reason why the parameter t\ does not appear 
in the final expression for the pair-production rate in a 
constant-crossed field is related to the fact that t\ is not 
a constant of motion (with respect to the classical equa¬ 
tions of motion, see Sec. IIIC). Correspondingly, the final 
(asymptotic) value of t\ will depend on the subsequent 
classical evolution of the charged particles and therefore 
on the nonlocal properties of the background field (see 
also the discussion in Sec. IllF). To avoid this problem, 
t-\ is not specified for a constant-crossed background field 
and only pair-production rates (not total probabilities) 
are calculated. For a laser pulse with finite duration, 
however, one obtains a total probability rather than a 
rate. Furthermore, as the particles leave the background 
field at a certain point, it is possible to specify the mo¬ 
mentum distribution of the final particles with respect to 
the parameter t\. 

Finally, we note that the amount nk M of absorbed laser 
four-momentum is related to the introduced Lorentz in¬ 
variant momentum parameters as follows 

1 I,, 2 

Pi +P 2 = ^ + nk^, n = -w — (1 +tl + t 2 ). (7) 

We will later show that this (asymptotic) value con¬ 
tains both a classical and a quantum contribution (see 
Sec. Ill E). As the laser-photon energy is only well defined 
for monochromatic fields, n is in general not an integer 
and only nk M is a meaningful quantity. Nevertheless, we 
will call n the number of absorbed laser photons. 


B. Total and differential probability 

Using the parameters introduced in Sec. II A, the total 
probability W (g, e) for the decay of a gamma photon with 
four-momentum q^ ( q 2 = 0 ) and polarization four-vector 
(eg = 0) into an electron-positron pair (see Fig. 1) can 
be written as (see, e.g., [ 8 ] and the literature cited in the 
introduction) 


W(q,e) 

d 3 W 

dRdt\dt 2 


, + 1/2 
/ dR 
J - 1/2 

TO 2 W 1 
( kq ) 2 8 (27t) 3 


f+°° . , d 3 W 

-oo ^ dRdtidt-2 ’ 


Y l-M(Pi>P 2 ;g )| 2 


spin 


( 8 a) 

( 8 b) 


where iM(pi,p 2 ;q) = ^nU Pl Q^(pi,q,-p 2 )v P2 is the re¬ 
duced iS-matrix element for the process [97] and we sum 
with respect to the final spin quantum numbers of the 
created particles (for simplicity they are not indicated). 
Here, u pi and v P2 denote the Dirac spinors for the elec¬ 
tron and the positron, respectively, and the nonsin¬ 
gular part of the laser-dressed vertex [see Eq. (10)]. 

In Eq. ( 8 ) the phase-space integrals are written in 
terms of the invariant momentum parameters defined in 


Eqs. (3) and (5) and d 3 W / (dRdtidt 2 ) represents the dif¬ 
ferential pair-creation probability with respect to those 
parameters. In order to calculate it we note that 

Y l- A/1 (Pi>I 5 2 ;g)| 2 = e M e* tr (y M (pi, g, —p 2 ) 

spin 

x{f 2 -m)g v (p 1 ,q,-p 2 )(f 1 +rn). (9) 

For on-shell momenta the nonsingular part of the dressed 
vertex is given by (see [ 8 , 92] for more details) 


G p (pi,q, -P 2 ) = (-*e)i ©0 g pp + Y ( G i®+i / 


[ip 

3 


3= 1,2 


+ G 2 © ii2 /| w )] +Hm7 5 Y f -'3©M/;"" • (10) 

j=l,2 ' 


where 


G 1 (R,kq)=e^, G 2 {R,kq) = ~ W 


kq 


2 (fcg) 2 


( 11 ) 


Furthermore, the so-called master integrals ©0 = 

&o(w,ti,t 2 ) and ©^ = &j i(w,ti,t 2 ) are given by (the 
notation agrees with [94]) 


r+ oo _ 

© 0 = / d(j> e iSr ^ w ' tl,t2 ’^\ 

( 12 a) 

r+ oo 


&j,i = / d 0 [^( 0 )] , e iSr < ,o > t 

J — OO 

( 12 b) 


with the nonlinear, field-dependent phase 

2 

— VO 777, 

Sr{w,tiM\4>) = S r (ti,t 2 ; 0 ), (13a) 

2 kq 


©r(U> Wi 0) — (1 + + t 2 )4> 

+ Y [ d $ 

i=l,2 J -°° 


(13b) 


From Eq. (10) we conclude that the pair-creation prob¬ 
ability can be calculated algebraically once the mas¬ 
ter integrals are known. As the integration range of 
<5j : i(w,ti,t 2 ) is naturally bounded, a numerical calcu¬ 
lation is readily accomplished (for l 7 ^ 0 , more details are 
provided in Sec. IIC). To determine © 0 {w,t\,t 2 ) we in¬ 
tegrate by parts and after neglecting the boundary terms 
we obtain the relation 


©o(w,fi,f 2 ) 


1 m 2 

-- - :—W 

2 n kq 


Y 2(w,t u t 2 ) 
2 = 1,2 


- 2 tj^j©j j i(ia, ti,t 2 ) 


(14) 


[n > 0, see Eq. (7)]. 
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C. Fourier-transformed master integrals 


From Eq. (13) we infer that the dependence of the mas¬ 
ter integrals &j t i(w,ti,t 2 ) on w is very simple. As a 
consequence, their Fourier transforms &jj(z,ti,t 2 ) [de¬ 
fined in Eq. (15)] can be calculated analytically. To this 
end we have to consider <5jj(w,ti,t 2 ) as a function of 
w £ (—oo,+oo), even though only the parameter range 
w £ [4, oo) is important from a physical point of view 
[the master integrals are everywhere well defined, see 
Eq. (12)]. After interchanging the order of integration, 
we obtain the following representation 


/• + 00 


®oA z MM)= / dwe l * w *=* 2 


= 47r 


m 2 ' 


(15) 


where the prime denotes the derivative with respect to 
the laser phase <f> [see Eq. (13)] and <j>~ is the (unique) so¬ 
lution of the equation = z. The unique¬ 

ness of <fi z follows from the fact that 


= 1 + X] [ l i ~ ( 16 ) 

i- 1,2 


is always greater than zero on the real axis. Thus, the 
calculation of the Fourier-transformed master integrals 
&j,i(z,ti,t 2 ) reduces to a root-finding problem (which is 
solvable numerically with low computational costs). 

Once <8j t i(z, ti,t 2 ) is calculated on a grid in z with suf¬ 
ficient resolution, the quantities which are 

related to &j.i(z,ti,t 2 ) by an inverse Fourier transform 
[see Eq. (15)], can be calculated numerically on a grid 
in w very efficiently by means of a single Fast-Fourier 
Transform (FFT) [98, 99]. Therefore, this approach re¬ 
duces the problem of calculating ti,t 2 ) as a func¬ 

tion of w, t\ and t 2 on a three-dimensional grid to an ef¬ 
fectively two-dimensional problem [from the viewpoint of 
computation costs, assuming that the root-finding prob¬ 
lem in Eq. (15) causes no significant overhead with re¬ 
spect to the FFT], In comparison with a direct calcu¬ 
lation of &j ; i(w,ti,t 2 ) [see Eq. (12)] using standard al¬ 
gorithms for highly, nonuniformly oscillating integrals, 
Eq. (15) reduces the required computational effort sub¬ 
stantially. 

Alternatively, one could also perform the change of 
variables </>—)• z = ©r(u>, ti, £2; </^) in Eq. (12) and evalu¬ 
ate the master integrals directly via FFT. The change of 
variables is one-to-one as & r (ti, t 2 \ </>) > 0, see Eq. (16). 
This approach has been applied in [100] to the analogous 
problem of nonlinear Thomson scattering. 


III. SEMICLASSICAL PICTURE 

By combining all relations presented in Sec. II, the nu¬ 
merical evaluation of the leading-order S-matrix element 


for the nonlinear Breit-Wheeler process and the determi¬ 
nation of the momentum distribution for the created par¬ 
ticles is straightforward. However, we obtain no further 
physical insights into the pair-production process in this 
way, as the S-matrix does not reveal any details about 
the dynamics taking place inside the interaction zone. 
Therefore, an intuitive semiclassical picture is now de¬ 
veloped, which is applicable for strong background fields 
(£ 1). Using optical lasers (photon energy lj ~ leV) 

intensity parameters £ > 100 are accessible at existing 
and upcoming laser facilities [9-13, 101]. In this regime 
the actual transformation from light to matter (which 
happens within a microscopically small formation region 
8<j) ~ l/£ in the laser phase) can be separated from the 
subsequent classical propagation of the created particles. 

In order to verify the reliability of the semiclassical 
approach (which is to a large extend similar to the one 
used in PIC codes), we compare its predictions with full 
numerical calculations. To this end we consider a linearly 
polarized laser field [ip = V’l, ^2 = 0, £ = £i] with the 
following pulse shape 

= sin 2 [0/(27V)] sin(0 + <£o) (17) 

for (j> £ [0, 27rlV] and zero otherwise. Here, N denotes the 
number of cycles and 4> 0 the CEP of the pulse (the nu¬ 
merical values of the parameters used in the calculations 
are given in the captions of the figures). Furthermore, 
we assume that the incoming photon has parallel polar¬ 
ization = A) 1 , see App. B). In this case the trace in 
Eq. (9) is given by 

-4Ai M A ll/ tr[...p 

e z 

= 2 m 2 (w - 4)[ - £ 2 |0i,i| 2 + 2£iiR(05<8 1 ,i)] 

+ 4m 2 |© 0 | 2 [2 1\ - (iu/2)(l +tl+ t 2 )] (18) 

and we denote the probability by W\\(q) = W(q, e = Ai) 
[see Eq. (8)]. 


A. Stationary-phase analysis 


To obtain an intuitive semiclassical picture, we apply 
now a stationary-phase analysis to the master integrals 
defined in Eq. (12). However, our calculation does not 
precisely follow the method of steepest descent, which 
is the standard approach (see, e.g., [37, 43, 49] and also 
[102]). Instead, it is shown that the integral along the 
real line is dominated by those points where the second 
derivative of the phase vanishes. Of course, the final 
result agrees with the one obtained using the method of 
steepest descent, but the derivation is less complicated 
as it does not require the deformation of the integration 
contour within the complex plane. 

From Eq. (13) we conclude that in the regime £ 3> 1 the 
master integrals are in general highly oscillating. From 
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the derivative of the phase 


2 

— W 171 

Sr(w,t 1 ,t 2 ]<j>) = 6rOi,t 2 ;0) (19) 

we infer that the master integrals have no stationary 
point on the real integration line (the prime denotes the 
derivative with respect to the laser phase <f>). Focusing 
on the case of linear polarization [see Eq. (16)] 

&r(ti,t 2 )<l>) = l+tl+[h -£^(<7 b )] 2 , (20) 

we find that the stationary points ip k of the phase 
Sr(<j>) = Sr(w, t\,t 2 ] <t>) [defined by = 0 ] are com¬ 

plex and given by 

i’i'Pk) = | {*1 ± V 1 + tij • (21) 

To obtain the leading-order approximation to the mas¬ 
ter integrals in the regime £ 1 , one could apply the 

method of steepest descent, i.e. deform the integration 
contour inside the complex plane such that it passes 
through the stationary points (see, e.g., [37, 43, 49]). 
However, the desired result is derived much faster by 
noting that the stationary points ip k are located pair¬ 
wise very close to the real line if < 1 (analogously 
to nonlinear Compton scattering and other processes 
[46, 49, 103] inside a plane-wave field; for |f 2 | 1 the 

pair-production probability is exponentially suppressed, 
see below). Mathematically, we have to deal with two 
stationary points which coalesce in the limit £ —> oo. For 
£ 1 the two stationary points nearly coalesce, a situa¬ 

tion which is discussed, e.g., in [104] (see also [94], App. H 
and [105], Chap. 36). 

Due to the presence of two stationary points ip k close 
to each (quasi-) stationary point (f> k defined by 

V#fc) = h/C, (22) 

we expect that the integral along the real line is domi¬ 
nated by small formation regions 5(f> around the points 
cj>k [for a linearly polarized laser field; in general we ob¬ 
tain two equations ipii't’k) — U/£i (* = 1 , 2 ) which should 
be fulfilled simultaneously]. From Eq. (20) we conclude 
that the oscillation frequency of the phase is as small 
as possible at these points (i.e. = 0 ) because the 

dominating contribution to the oscillation frequency be¬ 
comes stationary (an intuitive physical interpretation for 
this condition is given in Sec. IIIC). In the following, we 
will call the points <pk (and not ip k ) stationary points 
for simplicity (to stress the difference, the points <p k . are 
called true stationary points). Moreover, we assume that 
all stationary points <j>k are located sufficiently far away 
from each other, i.e. we ignore subtleties arising around 
the extremal points of [note that pair production is 
ineffective in these regions, as is small]. 

As the main contribution to the master integrals arises 
from the regions around the phases <pk where t\ ss 


£ip(4>k)i we expand the phase Sr{4>) = Sr(w,ti,t 2 ;<j)) 
[see Eq. (13)] around a stationary point 4>k up to cubic 
order 


Sr(<A) ~ S r (<j>k) + a(<j> - 0 fe ) + -b((f> - <f> k f 

2 


w m 
“= 2 t ^ (1 


4 ), 




(23) 


Here, we focus on the regime y > 1 (where the pair- 
production probability is not exponentially suppressed, 
see below) and define the formation region 5(f> = <fi — </>*, 
around </>*, by the requirement that the phase in Eq. (23) 
remains of order one (see, e.g., [106, 107]). As a ~ £ 
and b ~ £ 3 , the formation region scales as Sip ~ l/£ (we 
always assume £ 1 in this section). Correspondingly, 

both the linear and the cubic term have to be taken into 
account. Higher-order terms do not change the behavior 
of the phase significantly (within the formation region) 
and can be neglected to leading order. We will show in 
Sec. HIE that the scaling 5<j) ~ l/£ for the formation 
region is reasonable from a physical point of view, as 
the energy absorbed classically from the background field 
within 5<j> is sufficient to bring the particles on shell. 

After the change of variables from <j> to t = y/b{<j> — <j> k ) 
the phase is approximately given by [see Eq. (23)] 

Sv{4>) ~ Sr(4>k) +xt+ ^f 3 , (24) 


where 


' rb 

i 

Tb~ 


V 2 


_\x(<pk) I. 


2/3 


2 kq 


V 2 


(1 + t%), 

1 2/3 


I x(<t>k)\ 


(25) 


and the absolute value of %( 0 ) = xV’ / (^ > ) denotes the 
local value of the quantum-nonlinearity parameter \ = 
(kq/m 2 )^ [32, 37]. 

Finally, we obtain from the stationary point (f> k [defined 
by Eq. (22)] the following contribution to the master in¬ 
tegrals [in the case of a linearly polarized laser field, see 
Eq. (12)] 




kq 2 w/2 

m 2 w [\x(<t>k)\ 



2tt Ai (p) e iSr{4>k \ 

(26a) 


» / . ti , . . f kq 2 

« —®o (w,t 1 ,t 2 ) - l ( —2 

^ \m z 

27rAi , (p)V’ , (^fe)e iSr(?i '' ) , (26b) 


u/2 14/3 




where p = {w/[2 |x(^)|]} 2/3 (1 + t%), x{<t>) = X^'{4>), 
X = ( kq/m 2 )i .; and Ai denotes the Airy function [105]. 
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FIG. 2. Momentum distribution for the created electron-positron pair [see Eq. (8)] for the parameters y=l, £ = 5, A = 5 
and (f>o = 7r/2 [the longitudinal momentum characterized by R is integrated numerically and the incoming photon has parallel 
polarization (e^ = A^)]. The parameters £ = 5 and \ = 1 could be obtained by colliding 17 GeV photons head-on with 
optical (uj = 1.55 eV) laser pulses having an intensity of 10 20 w /cm 2 (note that few-cycle laser pulses are envisaged, e.g., at the 
petawatt field synthesizer (PFS) in Garching [108]). a) Full numerical calculation of the spectrum [see Eq. (15)]. b) Local 
constant-crossed field approximation applied on the amplitude level [see Eq. (26)] . The inset shows that the interference pattern 
is lost if the local constant-crossed field approximation is applied on the probability level, c) Outline for t 2 = 0. Solid line: full 
numerical calculation; dotted (dashed) line: local constant-crossed field approximation applied on the amplitude (probability) 
level. 


As the properties of the Airy function imply that pair 
production is exponentially suppressed for ^ < 1 and 
|t 2 1 ,w 1 [8, 32], we will consider % = 1 in the nu¬ 
merical calculations. Experimentally, the regime \ > lj 
£ 1 is accessible with presently available technology, 

i.e. by colliding GeV photons (obtainable, e.g., via Comp¬ 
ton backscattering [109-116]) with strong optical laser 
pulses. 


B. Interference substructure of the spectrum 

To compare the results derived in Sec. Ill A with those 
already known in the literature, we consider first a 
constant-crossed background field [^(</>) = (j)\. In this 
case the stationary-point equation t\ = £,i/j((f>k) = 

[see Eq. (22)] has only one solution and the approxima¬ 
tions leading to Eq. (26) are exact. Correspondingly, we 
obtain the probabilities (rates) for pair creation inside a 
constant-crossed background field given, e.g., in [37, 39] 
by combining Eq. (18) with Eq. (26). 

For an oscillatory plane-wave background field, how¬ 
ever, one finds in general more than one stationary point. 
Physically, this implies that the electron-positron pair 
can be created with the same asymptotic quantum num¬ 


bers at different laser phases (note that different forma¬ 
tion regions are usually separated on the macroscopic 
scale given by the laser wavelength). In accordance with 
general principles we expect that the existence of differ¬ 
ent pathways with the same final state causes interference 
effects similar to those in a multi-slit experiment (for the 
importance of interference effects see also [45-53]). 

The presence of interference fringes in the spectrum is 
demonstrated in Fig. 2. The first subplot was obtained 
using a full three-dimensional numerical calculation of 
the spectrum based on the method introduced in Sec. IIC 
[see Eq. (15)]. As shown in the second subplot, the semi- 
classical approximation introduced in Sec. Ill A is already 
for £ = 5 in very good agreement with the exact result. 

Due to the fact that the background field is approxi¬ 
mated locally around each stationary point as a constant- 
crossed field during the derivation of Eq. (26), we call the 
semiclassical approximation also local constant-crossed 
field approximation. However, it is important that the 
local constant-crossed field approximation is applied on 
the level of the probability amplitude (i.e. we do not 
simply average the spectrum of a constant-crossed field 
over the laser pulse shape). The essential difference be¬ 
tween both approaches is the presence of the phase factor 
exp [iSr{<pk)] hr Eq. (26), which gives rise to the interfer- 









ence substructure once the contribution of multiple sta¬ 
tionary points is taken into account. If the spectrum is 
calculated for each stationary point separately and the 
result is added on the probability level, the interference 
fringes are lost (see inset in Fig. 2; so far this approach 
was called local constant-crossed field approximation in 
the literature). 

As the interference pattern is determined by the phase 
factor exp [iSr(0k)]> we conclude from Eq. (13) that the 
oscillation frequency of the interference fringes in the 
spectrum scales as ~ £ 3 for w, ~ £ 2 for t\ and ~ £ 
for t- 2 . Here, we define the oscillation frequency of the 
spectrum with respect to a momentum parameter x as 
the inverse of the change Sx which is needed to advance 
from one local maximum of the differential probability to 
an adjacent one. The difference between the oscillation 
frequencies for t\ and £2 is clearly visible in Fig. 2. 

In order to fully resolve the interference substructure 
of the spectrum we used for £ ~ 10 a grid in momen¬ 
tum space (w, ti, t 2 ) with ~ 10 5 x 10 4 x 10 3 = 10 12 
data points (to obtain the two-dimensional plots we inte¬ 
grated numerically over the remaining momentum vari¬ 
able). From the above scaling laws for the oscillation 
frequency we conclude that this choice ensures enough 
sampling points per oscillation period of the interference 
fringes. As a cross-check we ensured that the total pair- 
creation probability calculated here by integrating nu¬ 
merically over the complete spectrum agrees with the 
one obtained in [8] using the optical theorem. 

As the interference substructure is an intrinsic nonlo¬ 
cal effect, it cannot be included easily into existing PIC 
schemes. However, the exact resolution of the transverse 
momentum components is beyond the achievable preci¬ 
sion of existing codes. Therefore, their overall precision 
should be increased first before interference effects can 
be studied. 


C. Classical interpretation of the stationary points 

To obtain an intuitive interpretation of the stationary 
points discussed in Sec. Ill A, we consider the classical 
equations of motion for an electron (positron) inside a 
plane-wave laser field. They predict that the time evo¬ 
lution of the electron four-momentum P M (<^>) is given by 
[32, 117-119] 

/>"(« = ps + A- r-(4>,<h)Po, 

+ 2< 27) 

where Pq = P M (<(>o) denotes the four-momentum at the 
laser phase 4>o and 

J (j) 0 


= frbPiit) - M&o)] (28) 

i- 1,2 

the integrated field tensor [compare with Eq. (1)]. To 
obtain the corresponding result for a positron we have 
to change the sign of the charge (e —> — e) in Eq. (27). 
Note that a laser field does not have a dc compo¬ 
nent [^(±00) = 0], which implies that the electron 
(positron) four-momentum does not change asymptoti¬ 
cally [P M (—00) = P M (+oo)]. This observation is in agree¬ 
ment with the Lawson-Woodward theorem [120, 121], 
which states that a plane-wave laser field cannot acceler¬ 
ate particles. 

The classical time evolution becomes particularly 
transparent if the four-momentum is expanded in the 
canonical light-cone basis [see Eq. (A4)] 

P M (</>) = p{<j>)q M + a(4>)k^ + to ^2 T i(<fr)K- ( 29 ) 

i=l,2 

The conservation of kP = kPo implies that also p(cf>) = 
p(4> 0 ) = kP/kq is conserved and 

1 777 ^ 

o{4>) = 0 kP~ Q t 1 + ^ + 72 ^ ( 30 ) 

is determined by the on-shell condition P 2 {(f >) = m 2 . 
Therefore, the nontrivial dynamic is entirely described 
by the transverse degrees of freedom 

n{4>) = n(<i> 0 ) - €i[ipi(<f>) - V’i(^o)]- (31) 

Correspondingly, also T 2 {(j)) = T2 (</>o) is conserved for a 
linearly polarized laser field = 0]. 

Using the above results we consider now the clas¬ 
sical evolution of the electron and the positron four- 
momentum Pi (</>) and 7*2(0 ) j respectively. Requiring the 
boundary conditions pf(00) = p,; [see Eq. (3)], we obtain 


Pi(</>) = rV + s'{cj))k^ + m ^ ti(</>)Af, 

(32a) 

2=1,2 


P 2 O/O = -r(f - s((j>)k M ti(<t>) Af, 

(32b) 

i=l,2 


where r' = r + 1, 


=U - 

(33) 


[note that U = tj(oo) as ipi(oo) = 0] and [compare with 
Eq. (4)] 

"j 777 ^ 

SW= 2 (34a) 

1 777 ^ 

S ' (0) = 2P ^ [1 + *iW+*2W ] - (34b) 

Furthermore, the following relation holds [compare with 
Eq. (7)] 

Pi (<t>) +P 2 (</>) = = n((j))k^, (35) 
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FIG. 3. Left side: Numerically calculated differential pair-production probability as a function of the transversal momentum 
parameters fa and fa (the longitudinal momentum characterized by R is integrated numerically). The incoming photon has 
parallel polarization (e M = A] 1 ), the laser pulse N = 2 cycles (such short pulses are envisaged, e.g., at the PFS in Garching 
[108]) and x = 1. We compare two different CEPs and two different intensities: a) <po = 0, £ = 10; b) (po = 0, £ = 5; c) 
(po = 7r, £ = 10; d) (po = 7r, £ = 5. The solid white lines confine the phase-space region where the pair can be produced at a 
phase <j> with \ip'((p)\ > 0.5 and the dashed white lines indicate the transverse momenta for which the pair can be produced at 
a local field peak. After integrating over fa and fa we obtain for the total pair-production probability W\\ = 0.09 % (£ = 10) 
and W\\ = 0.045% (£ = 5), up to this precision it is independent of the CEP. Right side: Plot of the laser pulse shape [solid 
line: ip'(<j>), dashed line: tp(<p)]. 


where 


[see Eqs. (33), (36) and (19)] 


2 

111 777 

«(0) = -^[l+i?(0)+il(0)]. (36) 


Assuming that a) the photon with four-momentum q M 
transforms within the short formation region 6<p ~ l/£ 
around a given laser phase (p c into an electron-positron 
pair [see Eq. (23) and the discussion below] and b) the 
charged particles subsequently obey the classical equa¬ 
tions of motion, we conclude from Eq. (35) that the 
four-momentum Q% = n(0 c )fc M must be absorbed “non- 
classically” from the background field during the pair- 
creation process itself (i.e. within the formation region). 
As the direct transformation of a real photon into a real 
electron-positron pair is kinematically forbidden, n(0 c ) 
is always greater than zero. Stated differently, the four- 
momentum QO is needed to bring the massive parti¬ 
cles on shell with the right initial conditions Pi(<p c ) and 
P 2 (0c) such that the subsequent classical propagation re¬ 
sults in the asymptotic four-momenta = p±( 00) and 
P 2 = P 2 { 00 )- 

To verify the correctness of this semiclassical picture, 
we demonstrate now that it is in perfect agreement with 
the results obtained in Sec. Ill A (i.e. it is valid in the 
regime £ 1). To this end we note the following relation 


n{cp) 


w m 2 
2 kq 


{1 + 

i=l,2 


= S' T {w,t ll t 2 \(p). 


(37) 


Correspondingly, the quantity n(<p), which is defined in 
Eq. (36) based on the classical equations of motion, cor¬ 
responds exactly to the oscillation frequency of the mas¬ 
ter integrals [jSp(to, ii, fa; 0), see Eq. (19)], which are ob¬ 
tained from the full quantum calculation of the pair- 
creation probability [see Eq. (12)]. 

As discussed in Sec. Ill A, the master integrals are 
dominated by small formation regions 5(p ~ l/£ around 
the (quasi-) stationary points 0*, defined by fa = £i0i(0fc) 
[i.e. fa(0fe) = 0; see Eq. (22) and Eq. (33)]. This result of 
the stationary-phase analysis has a very intuitive expla¬ 
nation within the semiclassical picture developed above. 
According to Eq. (37) the stationary points correspond 
exactly to those laser phases where n(0) is minimal. As 
the non-classical absorption of the four-momentum Q^ 
can be depicted as a tunneling process (see also [14]), it 
is natural to interpret n(0) as a measure of the effective 
tunneling distance. Therefore, the above intuitive picture 
also predicts that the pair-creation process happens pre¬ 
dominantly at those laser phases where n(0) is minimal, 
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in agreement with the full quantum calculation. 


D. Scaling laws for the spectrum 

From the discussion in the previous section we expect 
that many properties of the asymptotic momentum dis¬ 
tribution can be understood from the classical equations 
of motion [see Eq. (29)]. As an example, they predict 
that the spectrum extends up to in the variables ti in 
the general case of elliptical polarization. 

This supposition is confirmed by the stationary-phase 
analysis carried out in Sec. Ill A. For the special case of 
a linearly polarized laser field [ip 2 {P) = 0] we obtain 
|fi| < £, as for \t±\ > £ the stationary-point equation 
fi = £p{pk) has no solution [see Eq. (22)]. However, the 
differential probability with respect to £2 is now entirely 
determined by quantum effects during the pair-creation 
process itself, which implies |<21 1 [^2 is a constant of 

motion if p 2 (p) = 0]. The correctness of this bound can 
be seen explicitly from the argument p of the Airy func¬ 
tion [see Eq. (26) and the discussion below]. The observa¬ 
tion that pair production is exponentially suppressed for 
|^21 1 is in agreement with the fact that the true sta¬ 

tionary points ip k are located far away from the real axis 
in this case [see Eq. (21)]. Both scaling laws are verified 
numerically in Fig. 3. 

Moreover, Fig. 3 shows that all qualitative features of 
the spectrum in t\ except the interference substructure 
can be understood from the classical acceleration of the 
charged particles in the laser field. In particular, the po¬ 
sition and extension of the spectrum (solid lines) is pre¬ 
dicted well by the classical equations of motion. Further¬ 
more, the highest pair-creation probability is obtained for 
those momentum parameters which require that the pro¬ 
cess happens around a peak of the laser intensity (dashed 
lines). 

We point out that for short laser pulses the classi¬ 
cal acceleration has a preferred direction, which depends 
strongly on the CEP of the pulse. This can be seen from 
the plot of p(p) in Fig. 3, which determines the final mo¬ 
mentum of the particles [see Eq. (31)]. Correspondingly, 
the large CEP effects visible in the spectrum can be un¬ 
derstood from classical physics (for the nonlinear Breit- 
Wheeler process they were first reported in [22]). 

For a linearly polarized background field \ip2{(p) = 0] 
R and t 2 are constants of motion [see Eqs. (5) and (32)]. 
Therefore, the differential probability distribution with 
respect to R and t 2 is entirely determined by the pair- 
production process itself and remains invariant under the 
subsequent classical propagation of the particles. This is 
demonstrated in Fig. 4. After integrating over t±, the 
spectrum looks very similar to the one obtained for a 
constant-crossed field (after averaging over the laser pulse 
shape). In particular, we observe no substructure due to 
interference effects. For a fixed value of t\, however, the 
differential spectrum shows clear interference fringes (see 
inset in Fig. 4). 


Finally, we note that the scaling laws given above for 
the momentum variables imply that the available phase- 
space in R 1 ti and t 2 for the nonlinear Breit-Wheeler 
process is proportional to £. Correspondingly, the lin¬ 
ear increase of the total pair-production probability as a 
function of £ in the regime £ 3> 1 [8] is a pure kinematic 
effect. 


E. Classical vs. quantum absorption of 
laser four-momentum 

The results obtained in Sec. Ill C [in particular 
Eq. (36)] allow us to distinguish theoretically be¬ 
tween the four-momentum which is absorbed quantum- 
mechanically from the laser field during the pair-creation 
process itself (quantum absorption) 

= p*p(p k ) + P 2 (Pk) - q M (38a) 

and the four-momentum which is transferred classically 
from the laser to the charged particles during the accel¬ 
eration of the particles (classical absorption) 

ridk^ =p l p+pP l - [p^ipk) +p 2 (p k )}- (38b) 

Asymptotically, however, we observe only the sum of 
both processes [see Eq. (7)] 

n = n q + n c \. (38c) 

Here, Pi(pk) and p 2 {pk) \pi 2 (Pk) = m 2 } denote the ini¬ 
tial values for the classical propagation, which starts at 
the laser phase pk once the real pair is already produced 
[see Eq. (32) and Sec. Ill F]. Correspondingly, rep¬ 
resents the four-momentum which must be transfered 
non-classically to bring the particles on shell during the 
creation process (see Sec. Ill C) and n c \k ^ denotes the 
four-momentum which is transfered after the produc¬ 
tion until the pair leaves the laser field with asymptotic 
four-momenta p*p and p respectively. Note that the 
classical acceleration of the particles (i.e. the absorption 
of n c 1 laser photons) does not contradict the Lawson- 
Woodward theorem [120, 121], as the charged particles 
are created inside the laser pulse. 

In the case of linear polarization \p 2 {p) = 0] t 2 is a con¬ 
stant of motion and the stationary points are determined 
from ti = £p{pk) [see Eq. (22)]. Therefore, 

1 ?n 2 9 . 1 to 2 9 .__ 

Uq= 2 + Ucl = 2 W ~k^ tl ' ( 39 ) 

Using the scaling laws discussed in Sec. Ill D we conclude 
that the quantum and the classical absorption scale as 
n q ~ £/x and n c 1 ~ £ 3 /x, respectively. 

An intuitive understanding of the scaling law for n q is 
obtained by squaring Eq. (38), which shows that 

2 

n q > 2 (40) 
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FIG. 4. Differential probability with respect to the parame¬ 
ters t 2 and R for x = 1, ^ = 10, IV = 5 and <f> o = 1^12 [full 
numerical calculation, t\ is integrated numerically and the in¬ 
coming photon has parallel polarization (e M = A(“)]. The inset 
shows d 3 W\\/{dRdtidti) for ti = 0 [in arb. units, see Eq. (8)]. 
The pronounced interference pattern vanishes after the inte¬ 
gral in t\ is taken. Note that the (spin-summed) differential 
probability does not depend on the sign of R and 


tained for a constant-crossed field using the 5'-matrix ap¬ 
proach agrees with the initial momentum distribution in 
R and t 2 immediately after the particles are brought on 
shell. Therefore, this approach is reasonable for the two 
parameters R and t 2 . 

In contrast, the appropriate initial condition for fi(</>&) 
is far from clear, as t \(</>*,) is not conserved by the classical 
equations of motion. At first sight one could expect the 
existence of a distribution function which determines the 
initial condition for t\ (similar as for R and f 2 ). However, 
we showed in Sec. Ill C that the fixed value ti(^fc) = 0 
should be used. This condition defines the stationary 
points [for linear polarization, see Eq. (22)] and therefore 
leads to Eq. (26) which represent the probabilities inside 
a locally constant-crossed field. Correspondingly, the us¬ 
age of Eq. (26) to determine the probability distributions 
for R and < 2 requires the initial condition ti (</>&) = 0 for 
the classical propagation. 


IV. SUMMARY AND CONCLUSIONS 


This lower bound becomes an equality at the pair- 
production threshold, which is obtained for Pi(4>k) = 
(he. for w = 4 and f 2 = 0 in the case of linear 
polarization). 

Note that although we indicated n q as the number of 
“quantum-absorbed” laser photons (referring to the fact 
that the electron and the positron are virtual inside the 
formation region), the scaling law n q ~ £ (for \ > 1) 
can be understood by noting that the four-momentum 
which is transferable classically during the formation re¬ 
gion 8<f> ~ l/£ [see Eq. (23) and the discussion below] 
scales as ~ £fc M [see Eq. (29) and Eq. (36)]. 

In conclusion, the classical energy transfer during the 
propagation of the particles is much larger than the en¬ 
ergy transfer which takes place during the pair-creation 
process itself in the regime £ 3> 1 (n c i n q ). Corre¬ 
spondingly, a possible depletion of the background field 
(e.g., due to the development of a QED cascade) is dom¬ 
inated by the classical energy transfer from the laser to 
the pairs. 


F. Initial conditions for the classical propagation 

The semiclassical picture established in Sec. Ill C 
agrees with the basic principles used in PIC codes. In or¬ 
der to include the Breit-Wheeler process into a PIC code, 
the main difficulty is to determine the correct initial con¬ 
ditions for the classical propagation of the charged par¬ 
ticles. To this end the spectrum obtained for a constant- 
crossed field is usually employed (see Sec. Ill A). As R 
and t ~2 are constants of motion [see Eqs. (5) and (32)], 
the asymptotic momentum distribution in R and f 2 ob- 


In the present paper the momentum distribution for 
electron-positron pairs produced via the nonlinear Breit- 
Wheeler process has been investigated for short laser 
pulses. Using a newly developed numerical scheme (see 
Sec. IIC) we have calculated for the first time the spec¬ 
trum on a three-dimensional lattice which fully resolves 
the interference substructure even in the ultra-relativistic 
regime £ 1. Furthermore, we have investigated the lo¬ 

cal constant-crossed field approximation and showed that 
it reproduces the spectrum (including the interference 
fringes) for £ 1 if it is applied on the probability- 

amplitude level. Correspondingly, three effects determine 
the final momentum distribution in the regime £ 1: 

The production of an electron and a positron with phys¬ 
ical mass m inside a constant-crossed field, their sub¬ 
sequent classical acceleration by the laser field and the 
interference between all production channels which lead 
to the same asymptotic quantum numbers. Accordingly, 
we verified that the produced electron and positron be¬ 
have like classical particles after they have left the for¬ 
mation region (see Sec. Ill, in particular Sec. Ill C and 
Sec. HID) and that the substructure of the spectrum 
can be attributed to interferences between the contribu¬ 
tions of different formation regions similar to those in a 
multi-slit experiment (see Sec. Ill B). Furthermore, it is 
shown that one can distinguish between a classical and 
a quantum absorption of laser photons (see in particu¬ 
lar Sec. HIE). As the former is dominant in the regime 
£ 1, a possible depletion of the laser field during the 

development of a QED cascade is mainly caused by the 
classical acceleration of the created charged particles. In 
summary, the new findings presented here allow a disen¬ 
tanglement between classical and quantum aspects of the 
pair-production process. 
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Appendix A: Light-cone coordinates 


We call the four four-vectors k^, ef (i G 1,2) a 
light-cone basis if they obey [92] 


k 2 = k 2 = 0, kei = ka = 0, 
kk — 1, e^e^ — Sij. 


(Al) 


Using the above properties and the determinant identity 
for S'fS one g nc [ s an y suc h light-cone basis 

obeys fl 2 = 1, where 

n = e^A^efef (A2) 


As any set of four four-vectors ef (i € 1,2), k^ and 
which obeys the relations given in Eq. (Al) represents 
a light-cone basis, it is natural to ask which expressions 
are invariant under a change of the underlying light-cone 
basis. To this end we consider two different bases k tl , ef 
and k'^, e'- 1 and denote the corresponding components of 
a four-vector v M by 

v + = k^Vp, v 1 = ef'iv, ti 11 = e^Vfj, 

/+ Tin a /U /II /« 

= v =efv M , v =e£v IJj , 

(v~ = v'~ = kv). The three coordinates (— ,_L= i, n) 
define a closed subspace and we obtain the relation 



0 

-efei 

-e' 2 ei 


° \ 

-eie 2 

—e 2 e 2/ 



(A7) 


In order to show that its determinant has magnitude one, 
we write 


ef 1 = ae f + 6ef + A k^, 

, (A8) 

e(f = cef + def + /rfcU 

As ef 2 = ef 2 = — 1 and efef = 0, we obtain a 2 + b 2 = c 2 + 
d 2 — 1 and ac + bd = 0. Without restricting generality, 
we set a = cos p, b = siny? and d = cos 6, c = sind. 
Finally, we obtain the two solutions 6 = —ip and 6 = 
—<p + 7r, which correspond to ad — be = ±1. Therefore, 
the measure dv~dv ± = dv'~dv , ' L and the delta function 


is called the orientation of the basis. In light-cone coor¬ 
dinates the metric g^ = diag(+l, — 1, — 1, — 1) is given 
by 

g^ = k»k v + k^k u - efef - efef. (A3) 


Intrinsically, the photon decay inside a plane-wave field 
[see Sec. II, in particular Eq. (1)] is characterized by the 
two light-like four-vectors and q^ and the constant 
field tensors f? v . Therefore, it is natural to expand the 
four-momenta of the created electron-positron pair in the 
following light-cone basis 


g M , Af 


/rv 

kq 


Af 


kq a/^o! 


(A4) 


6^’ ± \v) = S(v-)S(v 1 )S(v u ) 

= S(v'~)5(v n )5(v ni ) = (A9) 

are invariant under a change of the light-cone basis. If 
v~ = d _l = 0 also the component v + = is invariant. 


Appendix B: Photon polarization density matrix 

The (complex) polarization four-vector e M of a photon 
with four-momentum q M ( q 2 = 0) must obey = 

— 1 and qe = 0. In the light-cone basis Af [see 

Eq. (A4)] we obtain for the metric [see Eq. (A5)] 


(we will always excluded the trivial case kq = 0 where 
pair production is forbidden and assume that kq ^ 0), 
which fulfills the completeness relation 

= JL (k»q v + q ^k v ) - Af Af - Af Af (A5) 

(the two four-vectors Af 2 have been used previously by 
several authors to analyze processes within strong plane- 
wave background fields, see, e.g., [129, 130]). To distin¬ 
guish the set of four four-vectors given in Eq. (A4) from 
other light-cone bases we will call it the canonical light- 
cone basis. 


g^ = -t- (fc V + q»k v ) - Af Af - Af Af 
kq 

and thus the polarization four-vector is given by 

= ciAf + c 2 Af + c. 3^, 

ke 

Cl = — (eAi), c 2 = — (eA 2 ), c 3 = —— 


(Bl) 


(B2) 


with the normalization condition |ci| 2 + |c 2 | 2 = 1. As 
the contraction of the matrix element with the four- 
momentum must vanish due to gauge symmetry, we 
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can restrict us to the vectors A/ and replace the density 
matrix by 


pM* = ^ £ pijA^Aj, 

i 4—1 2 

(B3a) 

1 i2 \ \2 


P 11 — Cl 1 P 22 — C 2 , 

Pl2 = C1C2, P21 = C[C2- 

(B3b) 


Correspondingly, it can be described by two Stokes angles 


si = cos ((p) sin(0), S2 = sin(y>) sin(0), 
s 3 = cos(0). 


(B7) 


Using the trigonometric identities 


The 2x2 density matrix p^j is Hermitian and has unit 
trace 


pij — A.if l kj v p'~ , plj — Pji — Pij , 

tr p= ^2 pa = 1. 

2 = 1,2 




J 


(B4) 


Any Hermitian 2x2 matrix can be expanded using the 
Pauli matrices & l and the identity 1 (with real parame¬ 
ters). Since tr a 1 = 0, we obtain (see [6], Eq. 8.9) 


P = 



+ SiCr *) 


1 / 1 + s 3 Si is 2 \ 

2 ySi + iS-2 1- S 3 J 

|ci| 2 Cic|A 

4^2 \C2\V ’ 



(B5) 


where Si are called Stokes parameters. The Stokes vector 
s = (si, S2, s 3 ) is a unit vector, which can be seen from 


^(1 + cos0) = cos 2 (0/2), 

^(1 — cos0) = sin 2 (0/2), (B8) 

2 cos(0/2) sin(0/2) = sin(0) 

we conclude that the complex coefficients ci and C2 can 
be expressed in terms of the stokes angles as 

Ci = cos(0/2) e~ %v ^ 2 , c 2 = sin(0/2) e +l¥,//2 , (B9) 

implying the representations 

|ci| 2 = cos 2 (0/2), |c 2 1 2 = sin 2 (0/2), 

1 (BIO) 

C1C2 = - sin(0) [cos(<p) — isin(<^)] 


det p = 0= -(l — 
4 v 


(B6) 


[note that we can always multiply by a total phase in Eq. 
(B9)]. 
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